THE GREAT ANNIHILATOR 1E1740.7-2942: MOLECULAR CLOUD CONNECTION AND 

CORONAL STRUCTURE. 



o 

Q 



> 

(N 

o 



Osmi Vilhu, Diana Hannikainen, Panu Muhli and Juhani Huovelin 

Observatory, University of Helsinki, Finland 

Juri Poutanen 

Stockholm Observatory, Saltsjobaden, Sweden 

Philippe Durouchoux 

DAPNIA, Service d'Astrophysique, CE Saclay, France 

Pierre Wallyn 

Jet Propulsion Laboratory, CalTech, USA 



ABSTRACT 



Using ^^CO and ^"^CO observations we present col- 
umn density maps of the molecular cloud (Ylsr = 
-135 km/s) in the direction of 1E1740.7-2942. Hydro- 
gen column densities of the cloud scatter between Nh 
= (3.5 - ll)xlO^^ cm~^, depending on the method 
used. From this we conclude, deriving first a simple 
analytic formula, that despite of the weakness of the 
iron fluorescent 6.4 keV line (Churazov et al. 1996) 
the source may lie inside the cloud, or at least close 
to its edge. The combined ASCA/BATSE spectrum 
from September 1993 and 1994 can be modelled with 
a two-phase accretion disc corona model, where the 
hot region is detached from the cold disc. Geomet- 
rically the hot phase can be interpreted e.g. as a 
number of active regions (magnetic loops) above the 
disc, or as a spherical hot cloud around the central 
object. 
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1. INTRODUCTION 



The discovery of a massive molecular cloud (at Vlsr 
« -130 km/s) in the direction of the bright X-ray 
source 1E1740. 7-2942 (Bally and Leventhal 1991; 
Mirabel et al. 1991) raised a picture of a compact 
source embedded in a molecular cloud. It was sug- 
gested that the unique properties of this object are 
related to the presence of dense gas surrounding it. 
In particular, the Bondi-Hoyle accretion onto a single 
black hole was discussed. 



The ASCA data imply N^ = 



low energy absorption (Sheth et al. 1996), and twice 
of that for the iron K-edge absorption (Churazov et 
al. 1996). Based on these column density estimates 
Churazov et al. concluded that the weakness of the 
iron 6.4 keV line suggests that the source is located 
behind the cloud. 



2. THE MOLECULAR LINE OBSERVATIONS 



The observations were performed between August 1 - 
8, 1993, using the Swedish ESO Submilhmetcr Tele- 
scope (SEST) situated at La Silla in Chile. The front- 
end was a Schottky receiver and the 1086 MHz band- 
width acousto-optical spectrometer (AOS) was used 
as the back-end. The intensity calibration was per- 
formed using the chopper-wheel method, the back- 
ground was eliminated with the position switching 
mode. The half power beamwidth (FWHM) of the 
telescope at 115 GHz is 45 arc sec. 

We observed two lines i2cO(J=l-0) (115 GHz) 
and i3C0(J=l-G) (110 GHz) at 35 spatial positions 
around 1E1740. 7-2942. In both transitions the hue at 
ViSfi = -135 km/s can be clearly identified (Fig. 1) 
and the antenna brightness temperature was inte- 
grated between —115 km/s and —150 km/s to give 
the total flux Fdoud- Three different methods were 
used to calculate Nh = 2xN(H2): 



Method 1. The standard conversion (see Bally and 
Leventhal 1991; Fig. 2) 

N(H2)/Fcioud(^^C0) = 2.6x10^0 cm'^/K km s-\ 

Method 2. Scaling with the LTE-relation given by 
Bally and Leventhal (1991) with T^xc = 15 K (see 
Fig. 3) 



N(H2)/Fcioud("C0) = S.lxlO^o cm-VK km s-\ 

Method 3. Assuming LTE like in the method 2, 

but computing Texc and r with the help of ^^CO 
and ^^CO, and finally using the conversion N(H2) = 
1.0xlO<5xN(i3cO) (Gahmetal. 1993; see Fig. 4). 

The column density maps based on these methods are 
shown in Figures 2-4. In the direction of 1E1740.7- 
2942, the methods 1, 2 and 3 give values of Nh = H, 
3.5 and 7, respectively (in units of 10^^ cm~^ and in 
the position of the circle centrum in Figs. 2-4). 
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Figure 1: CO (solid) and ^'^ CO (dotted) spectra at the 
(0,0) position of the maps in Figs. 2-4- 



3. IS 1E1740.7-2942 INSIDE THE CLOUD? 



Churazov, Gilfanov and Sunyaev (1996) showed that 
in the ASCA data there is no (or very weak) evi- 
dence o£ the fluorescent iron line at 6.4 keV. The line 
equivalent width is smaller than 20 eV with 90 per 

cent confidence, the 'best value' being around 10 eV. 
Further, Churazov et al. showed that the column 
must be less than (2 - 3)xl0^^ cm^^, if the source 
is in the center of the cloud. They concluded that 
1E1740.7-2942 is behind the cloud. 

Our measurements in the direction of 1E1740. 7-2942 

N/j = (3.5 -ll)xl0^^ cm~^ suggest that the source 
location may be behind the cloud to explain both 
the iron K-edge absorption and the standard IS- 
absorption in this general direction (5-7) x 10^^ cm~^ 
(Sheth et al. 1996). However, since the total col- 
umn (integrating over all velocities in Fig. 1) is large 
(6x10^3 cm-2), 1E1740.7-2942 may equally well lie 
in front of the cloud. 

Assume that 1E1740. 7-2942 is at a distance d from 
the centre of a spherical homogeneous cloud with ra- 
dius r. For a Thomson optically thin cloud, it is 
rather simple algebra to derive the following relation 
between the equivalent width EW of the fluorescent 
line, a = d/r, the column density across the cloud 
center Nh and iron abundance NFe (number den- 



Figure 2: The Hydrogen column density map around 
IEI74O.7-2942 from the ^'^CO line observations with a 
linear scaling (method 1, see the text). The ASCA res- 
olution circle (PSF diameter 2.9 arc min = 7 pc at 8.5 
kpc) centered at the radio- jets (Mirabel et al. 1992) is 
shown. The darkest region cojTcsponds to the largest col- 
umn density. The range of Nh is between (9 - 13)xl(f'^ 
cm~^ inside the circle. 



sity): 



EW = 30 Nh, 23 {NFe/NFeo) G{a) 



(1) 



where the equivalent width EW is in units of eV, 

column density is in imits of lO^'^cm^^, and the geo- 
metrical factor G(a) is given by the formula: 



G{a) 



1 H In 



2a 



1 



(2) 



To derive the numerical factor in equation 1 we used 
the same values for the absorption cross section at 
the iron K-cdge (3.5x10"^" cm^ per iron atom ) and 
for the fluorescent yield {W = 0.35) as Churazov et 
al. (1996) did. 

For a solar abundance of iron (NFe/NH = 10~^-^) , 
this equation predicts a = 1.1 (i.e. the source is close 
to the edge of the cloud) for the largest column Njj = 
11 X 10^2 cm-2 (method 1) and for the 'best' EW = 10 
eV derived by Churazov et al.. The smallest column 
from the method 2 (3.5x10^^) gives a = 0.1 (i.e. the 
source is close to the center of the cloud). Values of a 
range between 0.9 and 1.5 for NFe = 2NFeQ which 
can not be rejected due to the large observed iron 
K-edge. 




Figure 3: As in Fig. 2, but from the CO line observa- 
tions with a linear scaling (method 2, see the text). The 
darkest regions correspond to the largest column density. 
The range of Nh is between (3 - 4)'XiO^^ cm~^ inside 
the circle. 



Figure 4- As in Figs. 2 and 3, but from the LTE method 
using CO and CO (method 3, see the text). The dark- 
est regions correspond to the largest column density. The 
range of Nr is between (5 - 9)xl(P^ cm~^ inside the 
circle. 



The real column of the cloud might be somewhat dif- 
ferent from the values used. However, the size of the 
scattering induced area inside the ASCA resolution 
circle is small enough to make our estimates reason- 
ably good (see Figs. 2-4, extreme values of N/f 
range between (3 - 13)xl0^^ cm^^ inside the circles). 
After all, the largest uncertainties in column density 
estimates come from the method used. 



4. ASCA/BATSE SPECTRAL MODELLING 



Wo attempted spectral modelling using two-phase 
disc-corona models described by Poutancn and 
Svensson (1996). In these models, soft radiation from 
the accretion disc gets comptonized by hot electron 
(-positron) gas in the corona. Different coronal ge- 
ometries can be considered: slab, hemisphere, and 
cylinder. The effects of different viewing angles and 
separation of the hot active regions from the disc 
(covering factor g) can also be modelled. This is par- 
ticularly relevant to 1E1740. 7-2942, since the weak- 
ness of the iron line suggests that the reflection com- 
ponent is small. This means that the disc is viewed 
edge-on and/or the cold disc covers much less than 
27r solid angle as viewed from the X-ray source. 

The models with large g tend to have too small fluxes 
at BATSE energies. To compensate this one can de- 
crease the number of reprocessed soft photons enter- 



ing the corona (small (^-factor). By the energy bal- 
ance this increases the optical depth, allowing more 
scattering orders at BATSE energies. Geometrically 
small g means that the hot corona is detached from 
the disc (like apexes of magnetic coronal loops) or 
is situated in a central hole inside the classical disc 
(Shapho ct al. 1976, Narayan 1996). 

We used as input observations those from ASCA 
(Sheth et al. 1996) and from the simultaneous 
BATSE archive observations analyzed by the JPL 
earth occultation method (Wallyn et al. 1996). We 
do not attempt to give 'the best fit' model, just 
demonstrate one possible model. 

Figure 5 shows the average ASCA spectrum (Septem- 
ber 26, 1993 and September 9-12, 1994; Sheth et al. 
1996) fitted well with a power-law fit (a — 1.1 and 
Nij = 8 X 10^2 cm-2), together with the BATSE data 
averaged over 13 days around the ASCA observing 
dates. In addition, the BATSE standard state spec- 
trum is shown (high state minus low state in the 1989 
-91 data to subtract the apparent background con- 
tamination). The spectrum obtained in this way is 
almost identical to the SIGMA standard state spec- 
trum (Simyacv et al. 1991). The high energy rise 
can not be modelled with thermal models. It needs 
a separate physical component, and part of that may 
come from a possible diffuse radiation in the BATSE 
field of view. 

Models with hemisphere type active regions are over- 
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Ftgure 5: The average ASCA spectrum of 1E1740.7- 
2942 from September 26, 1993, and September 9-12, 1994 
(Sheth et al. 1996), together with the BATSE spectrum 
averaged over 13 days around these dates (data points 
with error bars). The standard state BATSE spectrum 
(almost identical to the standard SIGMA spectrum by 
Sunyaev et al. 1991) is also shown (stars). Edge-on 
models with no reflection component and with detached 
active regions are overplotted (g — 0.05 and 0.01). For 
comparison, the same model with active regions in touch 
with the disc is shown by the solid line (g — 0.7). The x- 
axis is logioe where e = E(keV)/511 and the y-axis eguals 
to logio{e^ Fpfi) where Fph is the observed photon flux 
(ph/cm^ /sec/keV). 



plotted in Figure 5. They have coronal temperature 
Tcor = 100 keV and soft photons black body tem- 
perature Tbb = 50 eV. The models are edge-on and 
the reflection component was switched off due to the 

weakness of the iron line. Different geometrical fac- 
tors were used (e.g. g = 0.01 corresponds to the case 
when 1% of the radiation reprocessed and reflected 
from the cold disc enters back the active region) . The 
optical depths for the models with g = 0.7, 0.05 and 
0.01 are 1.05, 2.8 and 3.75, respectively. 

In the case of electron-positron pair dominated 
corona, the compactness parameter l^iss = 
(L/ H){aT /n^eC^) can be used to estimate the dimen- 
sion H of the active region. Assuming the luminosity 
L = 3 X 10'^'' ergs/s for the standard state (Sunyaev 
et al. 1991), one obtains for the best model in Fig. 5 
with g = 0.01 (Uss = 1046) H = 8x10^ cm. This 
corresponds to 3R5ch for a solar mass central object. 
This type of central hot flow solutions were found 
already by Shapiro et al. (1976). For hotter Tcor 
this dimension increases to match better with a more 
massive BH. Assuming many active regions, H scales 
inversely with the number. For e+e~ -pairs in the ac- 
tive region, the number density can be estimated in 
the model as n = 1.5 x W'^^t/H w 7 x lO^^cm-^. 
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